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ABSTRACT 


Weapon  systems  are  comprised  of  parts  that  are  subject  to  random  failures.  When  a  part 
fails,  it  must  be  replaced  by  an  operable  part  that  is  provided  by  a  supply  network  that 
supports  the  system.  Supply  networks  consist  of  many  locations  where  spare  parts  are 
held,  known  as  “echelons.”  Examples  include  depots,  fulfillment  centers,  and  customers. 
When  many  identical  weapon  systems  operate  in  parallel  and  rely  on  a  multi-echelon 
supply  network  for  replacement  parts,  decision  makers  must  choose  where  and  how  to 
invest  their  resources  into  the  purchase  of  spare  parts.  This  thesis  uses  stochastic 
optimization  to  leverage  those  decisions  in  order  to  maximize  the  expected  number  of 
time  periods  a  set  of  weapon  systems  are  available  for  use.  Specifically,  we  develop  a 
model  that  determines  the  optimal  stock  levels  of  spare  parts  to  store  at  each  echelon  of 
the  supply  network.  The  formulation  integrates  part  failure  uncertainty,  transit  times,  and 
monetary  constraints.  Model  outputs  also  provide  decision  makers  with  a  clear  estimate 
of  marginal  availability  gains  for  each  dollar  invested  in  purchasing  spare  parts. 
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EXECUTIVE  SUMMARY 


Weapon  systems  are  comprised  of  parts  that  are  subject  to  random  failures.  When  a  part 
fails,  it  must  be  replaced  by  an  operable  part  that  is  provided  by  a  supply  network  that 
supports  the  system.  Supply  networks  consist  of  many  locations  where  spare  parts  are 
held,  known  as  “echelons.”  Examples  include  depots,  fulfillment  centers,  and  customers. 
If  the  supply  network  does  not  provide  replacement  parts  in  a  timely  manner,  weapon 
system  availability  suffers.  When  many  identical  weapon  systems  operate  in  parallel  and 
rely  on  a  single  supply  network  for  replacement  parts,  decision  makers  must  choose  stock 
levels — the  most  efficient  quantities  of  spare  parts  to  stock  at  each  echelon  in  a  supply 
network,  given  a  limited  budget,  so  that  system  availability  is  maximized. 

This  thesis  develops  a  stochastic  optimization  model  that  prescribes  optimal 
investments  in  spare  parts  for  a  weapon  system  throughout  a  multi-echelon  supply 
network.  This  model  also  provides  decision  makers  with  an  estimate  of  expected  system 
availability  given  such  investments.  Our  formulation  models  the  dynamic  interactions  of 
random  part  failures,  transit  times,  and  monetary  constraints,  all  of  which  impact  the 
supply  network’s  ability  to  maintain  weapon  system  availability. 

Our  model  confronts  the  above  problem  in  two  steps.  First,  we  incorporate 
uncertainty  into  our  formulation  by  randomly  generating  scenarios  of  time -based  part 
failures,  which  become  input  parameters  for  our  stock  level  calculations.  Second,  we 
model  the  evolution  of  our  inventory  levels  over  time  and  incur  a  penalty  in  our  objective 
function  if  a  system  is  unavailable  (i.e.,  if  the  supply  network  does  not  provide 
replacement  parts  in  a  timely  manner  for  that  system).  Model  outputs  provide  decision 
makers  with  optimal  stock  levels  and  a  clear  estimate  of  marginal  availability  gains  they 
can  expect  for  each  dollar  invested  in  purchasing  spare  parts. 

We  implement  our  formulation  and  show  that  as  supply  network  topology 
becomes  larger  and  more  complex,  more  investments  in  spare  parts  are  needed  to  ensure 
high  rates  of  weapon  system  availability.  Also,  we  find  that  marginal  availability 
increases  faster  with  additional  investments  in  retail  budget  than  in  wholesale  budget. 


Moreover,  we  show  that  in  some  cases,  we  may  obtain  similar  availability  for  different 
investment  levels.  For  our  specific  example,  a  decision  maker  can  expect  96  percent 
availability  by  optimally  spending  any  budget  between  $80  and  $95  on  the  procurement 
of  spare  parts  for  the  system. 

Additionally,  we  examine  the  variation  in  our  model’s  output  due  to  uncertainty. 
We  show  that  increasing  the  number  of  part  failure  scenarios  over  which  we  optimize 
increases  both  solution  quality  and  the  amount  of  time  required  to  find  an  optimal 
solution.  We  conclude  that  it  is  important  to  balance  both  aspects.  For  our  particular 
example,  we  illustrate  that  modeling  25  part-failure  scenarios  per  replication  gives  us 
results  that  reasonably  balance  objective  value  precision  and  solve  time.  Furthermore,  we 
show  that  as  the  number  of  part-failure  scenarios  modeled  increases,  our  model 
converges  towards  an  optimal  solution  for  our  wholesale  and  retail  stock  levels. 

Lastly,  we  examine  restriction  and  relaxation  techniques  that  decrease  the  time 
required  to  find  an  optimal  solution.  We  achieve  moderate  success  with  these  techniques 
but  conclude  that  further  research  in  this  area  is  needed. 
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I.  INTRODUCTION 


A.  BACKGROUND 

1.  Problem  Definition 

Weapon  systems  are  comprised  of  parts  that  are  subject  to  random  failures.  When 
a  part  fails,  it  must  be  replaced  by  an  operable  part  that  is  provided  by  a  supply  network 
that  supports  the  system.  If  replacement  parts  are  not  made  available  in  a  timely  manner 
by  the  supply  network,  the  weapon  system  becomes  unavailable  until  replacements 
arrive. 

Supply  networks  consist  of  many  locations  where  spare  parts  are  held,  known  as 
“echelons.”  Examples  of  echelons  include  depots,  fulfillment  centers,  and  customers. 
When  many  identical  weapon  systems  operate  in  parallel  and  rely  on  a  single  supply 
network  for  replacement  parts,  decision  makers  must  choose  where  and  how  to  invest 
their  resources  into  the  purchase  of  spare  parts  in  order  to  maximize  the  amount  of  time 
the  systems  are  available  for  use.  In  particular,  they  must  choose  stock  levels — the  most 
efficient  quantities  of  spare  parts  to  stock  at  each  echelon.  Should  decision  makers  choose 
stock  levels  that  are  too  low,  availability  of  the  weapon  systems  suffers  and  unnecessary 
strain  is  placed  on  the  supply  network.  Conversely,  should  decision  makers  choose  stock 
levels  that  are  too  high,  funds  are  wasted  or  inefficiently  used  (i.e.,  tied  to  parts  in  stock 
when  they  could  be  used  for  other  purposes).  As  a  result,  decision  makers  face  the 
problem  of  determining  how  to  optimally  choose  stock  levels.  In  particular,  given  a 
limited  budget,  a  decision  maker  must  determine  how  many  spares  of  each  part  to  buy, 
where  in  the  supply  network  to  stock  these  parts,  and  how  much  weapon  system 
availability  to  expect. 

To  address  these  questions,  this  thesis  formulates  a  mathematical  model  that 
determines  the  optimal  way  to  invest  limited  budgets  into  the  procurement  of  spare  parts 
to  maximize  the  expected  availability  of  a  weapon  system.  Specifically,  we  use  stochastic 
optimization  to  maximize  the  expected  number  of  time  periods  a  weapon  system  is 
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available  for  use.  The  key  decisions  in  our  model  determine  the  optimal  stock  levels  of 
spare  parts  at  each  echelon  in  the  system’s  supply  network. 

Our  approach  differs  from  other  research  in  that  most  supply  network 
optimization  studies  ignore  interactions  among  part  failures,  repair,  transit  times,  and 
costs.  Consequently,  decisions  made  to  improve  one  component  of  the  supply  network 
may  adversely  impact  other  components.  Without  taking  a  comprehensive  view  of  the 
system  and  making  decisions  accordingly,  one  cannot  hope  to  achieve  full  efficiency. 

The  key  contribution  of  our  research  is  the  integrated  modeling  of  uncertainty, 
time,  and  monetary  constraints  on  a  multi-echelon  supply  network.  We  incorporate 
uncertainty  into  our  modeling  by  randomly  generating  scenarios  of  time-based  part 
failures  for  a  particular  weapon  system.  Then,  we  use  those  scenarios  as  input  parameters 
for  our  stock  level  calculations  while  enforcing  non-anticipativity.  We  also  model  the 
evolution  of  our  inventory  levels  over  time  and  incur  a  penalty  in  our  objective  function 
if  the  supply  network  does  not  provide  replacement  parts  in  a  timely  manner.  Finally,  we 
provide  decision  makers  with  an  estimate  of  the  weapon  system  availability  they  can 
expect  given  an  optimal  investment  in  spare  parts  for  the  system.  To  our  knowledge,  no 
other  research  incorporates  these  three  aspects  into  a  single  mathematical  model.  This 
level  of  detail  comes  at  the  expense  of  lengthy  solve  times.  Even  for  very  small  supply 
networks  with  few  parts,  the  optimization  solver  requires  a  substantial  amount  of  time  to 
find  an  optimal  solution.  We  view  our  research  as  the  first  step  towards  high-fidelity, 
multi-echelon  modeling;  perhaps  future  research  can  build  upon  our  efforts  and  decrease 
the  computational  time  required  to  solve  our  models. 

2.  Weapon  System 

Figure  1  illustrates  a  notional  weapon  system.  The  system  is  comprised  of 
n  different  types  of  parts  pi,  p2,...,Pn-  These  parts  operate  in  series,  that  is,  the  system  is 
only  available  for  use  during  a  particular  time  period  if  all  parts  are  working  during  that 
period.  Each  part  has  a  known  probability  of  failing  after  operating  for  a  given  number  of 
periods.  We  assume  that  failures  for  each  part  are  independent  events  (i.e.,  failure  of  one 
part  does  not  impact  the  probability  of  failure  of  other  parts). 
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Figure  1.  Notional  weapon  system. 
a.  Weapon  System  Availability 

A  weapon  system  is  only  valuable  when  it  is  available  for  operational  use.  This 
occurs  when  the  quantity  of  working  parts  installed  and/or  spares  on  hand  is  at  least  the 
quantity  required  to  operate  the  system.  Accordingly,  our  research  maximizes  the 
availability  of  a  weapon  system  over  a  defined  time  horizon  by  finding  optimal  stock 
levels  of  each  part  for  each  echelon  in  a  system’s  supply  network.  We  define  availability 
of  a  weapon  system  as  the  fraction  of  time  that  system  is  available  for  use  and  calculate  it 
as  follows: 


availability  = 


^available 


^available  ^ not  available 


(1) 


where  and  ^..auaue  the  number  of  periods  the  system  is  and  is  not  available 

for  use,  respectively.  By  construction,  the  availability  metric  always  produces  a  value 
between  0  and  1,  providing  us  with  an  unambiguous  measure  with  which  to  compare  the 
effectiveness  of  different  solutions. 
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b.  System  of  Systems  Concept 


We  further  expand  the  scope  of  our  availability  metric  to  embody  the  “system  of 
systems”  concept,  illustrated  in  Figure  2.  We  reason  that  many  independent  and  identical 
weapon  systems  operate  in  parallel  in  one  large,  aggregate  system.  Accordingly,  instead 
of  measuring  the  availability  of  each  individual  system,  our  model  makes  stock  level 
decisions  that  maximize  availability  of  the  aggregate  system.  This  is  defined  as  the 
average  availability  of  the  individual  weapon  systems  that  comprise  the  aggregate 
system.  The  system  of  systems  concept  applies  to  our  problem  in  that  each  individual 
weapon  system  requires  replacement  parts  from  a  single  supply  network  that  stocks  spare 
parts  for  all  weapon  systems.  When  two  weapon  systems  require  the  same  part  at  the 
same  time  and  the  supply  network  only  has  one  spare  on  hand,  our  model  chooses  to 
issue  the  spare  to  the  entity  that  results  in  the  largest  average  availability.  If  we 
considered  each  weapon  system  in  isolation,  we  would  neglect  these  types  of  interactions 
and  improperly  estimate  the  aggregate  system’s  availability. 


Figure  2.  System  of  systems  concept. 


Mapping  Investments  in  Spare  Parts  to  Availability 


c. 


A  key  benefit  of  our  model  is  that  it  informs  decision  makers  of  how  much 
weapon  system  availability  they  can  expect,  given  monetary  investments  in  spare  parts 
for  that  system.  Figure  3  illustrates  a  notional  investment  versus  availability  curve.  The 
horizontal  axis  corresponds  to  our  total  monetary  investments  in  spare  parts  for  a 
particular  weapon  system  and  the  vertical  axis  corresponds  to  the  estimated  fraction  of 
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time  a  weapon  system  is  available  for  use.  Each  point  on  the  curve  is  a  notional  output  of 
our  model  corresponding  to  the  maximum  average  availability  value  that  can  be  achieved 
by  making  optimal  stock  level  decisions  given  a  particular  level  of  investment.  For 
example,  should  a  decision  maker  choose  to  invest  $20  into  the  purchase  of  spare  parts 
for  the  modeled  weapon  system,  the  curve  shows  that  the  weapon  system  is  expected  to 
be  available  60  percent  of  the  time.  Note  that  initial  investments  in  inventory  result  in 
large  availability  marginal  gains.  However,  as  monetary  investments  in  spare  parts 
increase,  these  gains  decrease  as  diminishing  returns  take  effect. 
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Figure  3.  Notional  curve  of  investment  in  spare  parts  versus  availability. 


3.  Repairable  Parts 

In  general,  parts  are  categorized  as  either  consumable  or  repairable.  Typically, 
consumable  parts  are  simple,  low-cost  items  that  are  discarded  after  failure.  Conversely, 
repairable  parts  are  expensive,  complex  items  that  may  be  discarded  after  failure,  but  are 
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often  repaired  by  a  servicing  echelon  and  placed  back  into  circulation  within  a  supply 
network.  Entities  that  operate  weapon  systems  comprised  of  repairable  parts  prefer 
purchasing  repaired  ones,  because  their  price  is  often  a  fraction  of  what  it  would  have 
cost  to  purchase  the  parts  as  new.  A  drawback  to  such  a  heavy  reliance  on  repairable 
parts  is  that  their  availability  within  a  supply  network  is  often  scarce  due  to  long  repair 
times.  Since  most  weapon  systems  are  comprised  of  repairable  parts,  our  analysis  focuses 
solely  on  these  parts. 

4.  Multi-Echelon  Supply  Network 

Figure  4  displays  a  simplified  multi-echelon  supply  network  for  a  repairable  part 
p.  Starting  at  the  bottom  of  Figure  4,  the  first  echelon  in  this  network  is  referred  to  as  the 
customer.  The  customer  installs  working  parts  to  operate  its  weapon  system  and 
maintains  a  stock  of  ready-for-issue  spares  to  replace  installed  parts  that  fail.  Moving 
upstream  from  the  customer,  the  second  echelon  is  referred  to  as  the  fulfillment  center. 
This  echelon  is  located  within  the  same  geographic  region  as  the  customer,  and  it 
maintains  a  stock  of  spares  that  are  used  to  supply  replacement  parts  to  the  customer. 
Finally,  moving  upstream  from  the  fulfillment  center,  the  third  echelon  is  referred  to  as 
the  depot.  The  depot  receives  failed  parts  from  the  customer,  repairs  them,  stores  them, 
and  issues  them  as  replacement  parts  to  supply  the  fulfillment  center.  All  events  take 
place  within  discrete  time  periods  we  denote  as  t. 
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Figure  4.  Basic  multi-echelon  supply  network  for  repairable  parts. 


a.  Multi-Echelon  Supply  Network  Events 

When  a  failure  of  part  p  occurs  at  customer  echelon  e  during  time  period  t,  the 
following  events  take  place: 

•  The  customer  checks  to  see  if  a  spare  p  is  on  hand: 

•  If  so,  the  failed  part  is  immediately  replaced  and  the 
weapon  system  remains  available  for  that  time  period.  The 
customer’s  quantity  of  on-hand  spares  is  decremented 
accordingly. 

•  Otherwise,  the  weapon  system  becomes  unavailable  for  the 
number  of  time  periods  that  elapse  until  a  replacement  is 
received  from  the  supporting  fulfillment  center. 

•  In  either  case,  the  customer  ships  the  failed  part  to  the 
supporting  depot  for  repair. 

•  If  the  customer  requires  a  replacement  part: 
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•  If  the  supporting  fulfillment  center  has  a  spare  on  hand,  it 
supplies  the  part  to  the  customer  and  its  quantity  of  on-hand 
spares  is  decremented  accordingly.  The  time  interval 
between  the  customer’s  part  failure  and  the  receipt  of  a 
replacement  part  from  the  supporting  fulfillment  center  is 
denoted  as  transit_time. 

•  Otherwise,  the  fulfillment  center  must  await  a  replacement 
part  from  its  supporting  depot.  The  number  of  time  periods 
it  takes  to  transport  a  replacement  part  from  the  depot  to  the 
fulfillment  center  is  also  denoted  as  transit_time. 

•  When  the  depot  receives  a  failed  part  from  the  customer,  it  repairs 
the  part  and  makes  it  available  for  issue  to  any  fulfillment  center 
that  requires  the  part.  The  number  of  time  periods  required  to 
repair  the  part  is  denoted  as  repair_time. 

b.  Enhanced  Multi-Echelon  Supply  Network 

Adding  to  Figure  4,  we  enhance  the  scope  of  our  multi-echelon  supply  network  to 
take  the  shape  of  Figure  5.  In  this  network,  multiple  depots  supply  replacement  parts  to 
multiple  fulfillment  centers,  which,  in  turn,  supply  replacement  parts  to  many  customers. 
We  assume  that  each  depot  specializes  in  repairing  particular  parts.  Consequently,  not  all 
depots  can  repair  all  parts.  Fulfillment  centers,  on  the  other  hand,  can  receive,  store,  and 
issue  all  parts.  Lateral  supply  only  exists  between  fulfillment  centers,  and  each  customer 
is  assumed  to  be  supported  by  only  one  (pre-specified)  fulfillment  center. 
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Figure  5.  Enhanced  multi-echelon  supply  network  for  repairable  parts. 


5.  Budget 

If  our  budget  were  unconstrained  and  we  had  unlimited  storage  space,  we  would 
choose  to  purchase  infinite  quantities  of  spare  parts  to  store  at  each  customer  echelon. 
This  would  ensure  that  aggregate  weapon  system  availability  would  always  be  100 
percent.  Unfortunately,  this  is  unrealistic.  Therefore,  we  must  constrain  our  model  by 
restricting  our  monetary  investments  in  spare  parts  (We  ignore  echelon  capacity,  which 
can  easily  be  incorporated,  if  needed).  Specifically,  we  include  two  constraints  in  our 
model,  reflecting  limits  on  our  expenditures  for  the  wholesale  and  retail  echelons  of  the 
supply  network.  The  wholesale  budget  represents  the  maximum  amount  of  funding 
available  to  procure  the  initial  quantity  of  spare  parts  to  be  stocked  at  depots  and 
fulfillment  center  echelons.  The  retail  budget  represents  the  maximum  amount  of  funding 
available  to  procure  the  initial  quantity  of  spare  parts  to  be  stocked  at  all  customer 


9 


echelons.  We  choose  to  model  wholesale  and  retail  budget  constraints  instead  of  a  single, 
network-wide  budget  constraint  so  that  our  formulation  aligns  with  real-world  budgeting 
practices  of  weapon  system  support. 

B.  LITERATURE  REVIEW 

To  establish  a  theoretical  framework  for  our  research,  we  build  upon  literature  in 
the  field  of  multi-echelon  modeling.  Work  in  this  area  dates  back  to  the  1950s  and  has 
grown  quite  extensively  in  recent  years  due  to  advances  in  computing  power.  We 
categorize  multi-echelon  research  into  deterministic,  stochastic,  and  simulation  models, 
and  borrow  key  ideas  to  guide  our  efforts. 

I.  Deterministic 

We  begin  by  examining  research  with  deterministic  parameters.  Sherbrooke,  often 
considered  a  pioneer  in  multi-echelon  network  modeling,  formulates  a  nonlinear  model 
that  calculates  stock  levels  by  minimizing  total  backorders  across  all  customers  [1].  He 
uses  a  technique  known  as  “marginal  analysis”  to  arrive  at  optimal  solutions.  Muckstadt 
expands  upon  Sherbrooke’s  formulation  and  uses  Lagrangian  relaxation  to  minimize 
expected  backorders  subject  to  a  constraint  on  inventory  investment  [2].  From  both 
authors  we  gain  an  understanding  of  the  link  between  weapon  system  availability  and 
monetary  investments  in  spare  parts.  Chandra  develops  a  dynamic  distribution  model 
with  warehouse  and  customer  replenishment  requirements  [3].  His  efforts  focus  on 
minimizing  transportation,  holding,  delivery,  and  order  costs  in  a  network.  From  his 
research,  we  gain  insight  on  how  to  formulate  inventory  balance  constraints  over  a  finite 
planning  horizon  and  how  to  model  the  evolution  of  a  supply  network  over  time.  Lee 
develops  a  continuous  review,  multi-echelon  model  for  repairable  items  when  emergency 
lateral  transshipments  are  allowed  between  customer  echelons  [4] .  His  efforts  expand  on 
lateral  shipment  work  previously  done  by  Sherbrooke.  Finally,  Pirkul  and  Jayaraman  use 
mixed-integer  programming  to  minimize  total  transportation,  distribution,  and  fixed  costs 
for  opening  and  operating  plants  and  warehouses  in  a  multi-commodity,  tri-echelon 
network  [5].  Customer  locations  and  their  demand  for  products  are  assumed  to  be  known 
in  advance. 
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2. 


Stochastic 


To  incorporate  uncertainty  into  our  modeling,  we  review  research  that 
incorporates  stochastic  processes.  Schneider  and  Kirkpatrick  introduce  us  to  stochastic 
optimization  and  show  us  approaches  for  heuristic  development  [6] .  Narrowing  our  focus 
to  supply-specific  research,  we  borrow  ideas  from  Axaster  on  how  to  model  lateral 
supply-determining  when  it  is  better  for  one  or  more  supporting  echelons  to  provide 
replacement  parts  to  a  supported  echelon  [7].  He  develops  decision  rules  for  lateral 
supply  in  a  single-echelon  inventory  system  consisting  of  a  number  of  parallel  echelons 
with  stochastic  demand.  His  approach  is  to  minimize  expected  costs.  Graves  determines 
inventory  stock  levels  in  a  multi-echelon  inventory  system  for  repairable  items  [8].  His 
work  presents  an  exact  model  for  finding  the  steady-state  distribution  of  net  inventory 
levels  for  parts  at  each  site  as  well  as  the  distribution  for  the  number  of  outstanding 
orders  for  a  site  at  a  particular  time.  From  his  work,  we  capture  the  events  associated  with 
modeling  repairable  parts. 

Tsiakis  et  al.  formulate  a  mixed-integer  program  that  determines  the  number, 
location,  and  capacity  of  warehouses  and  distribution  centers  to  be  set  up,  the 
transportation  links  that  need  to  be  established  in  the  network,  and  the  flows  and 
production  rates  of  material,  subject  to  uncertain  demands  [9].  They  seek  to  minimize 
total  annualized  costs  over  the  network.  From  their  work,  we  gain  an  appreciation  for 
approaches  on  how  to  handle  uncertainty  including  the  use  of  stochastically-generated 
scenarios.  Moreover,  we  align  our  objective  function  with  theirs,  in  that  we  want  to  find 
robust  solutions  that  perform  well  over  all  scenarios  of  uncertain  demand,  not  for  a  single 
presumed  outcome.  Caggiano  et  al.  describe  and  validate  a  practical  method  for  efficient 
computation  of  time-based  service  level  requirements  in  a  multi-item,  multi-echelon 
service  parts  distribution  system  [10].  They  use  channel  fill  rates  to  determine  where  to 
position  inventory  in  order  to  satisfy  customer  service  agreements.  Ganeshan  formulates 
a  multiple-retailer,  single-warehouse,  multiple- supplier  model  that  encompasses  both 
inventory  and  transportation  factors  [11].  Furthermore,  his  research  finds  near-optimal 
stocking  policies  for  reorder  points  and  order  quantities  under  stochastic  demands  and 
lead-time  constraints,  subject  to  customer  service  constraints.  From  his  work,  we  gather 
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ideas  on  the  interactions  between  inventory  and  transportation  decisions.  Neale  and 
Willems  provide  ideas  on  how  to  accommodate  situations  where  the  rate  of  demand 
changes  over  time  [12].  Their  formulation  determines  inventory  locations  and  levels  in  a 
supply  chain  facing  stochastic,  non- stationary  demand.  lida  considers  a  periodic-review, 
dynamic,  multi-echelon  inventory  problem,  with  non- stationary  demands  [13].  Moreover, 
his  work  expands  our  knowledge  base  by  showing  us  the  impacts  of  non-stationary 
demands  over  time.  Simchi-Levi  and  Zhao  provide  methods  for  evaluating  stochastic, 
multi-echelon  inventory  systems,  specifically,  the  queuing  inventory  method,  the  lead- 
time-demand  method,  and  the  flow-unit  method  [14].  The  queuing- inventory  and  lead- 
time-demand  methods  advise  us  of  system  performance  measures  at  random  points  in 
time,  whereas  our  approach  more  closely  resembles  the  flow-unit  method  in  that  it 
models  the  impacts  of  decisions  on  time -based  events.  Ettl  et  al.  formulate  a  constrained 
nonlinear  optimization  problem  that  minimizes  the  total  average  dollar  value  of  inventory 
in  a  supply  network,  subject  to  meeting  the  service-level  requirements  of  customers  [15]. 
Finally,  Acimovic  and  Graves  provide  us  with  ideas  on  how  to  enhance  the  functionality 
of  our  model  for  future  research.  Specifically,  they  formulate  a  model  that  minimizes 
average  outbound  shipping  costs  in  a  two-echelon  network  where  business  rules 
determine  what  fulfillment  centers  service  which  customers  [16].  From  their  work,  we 
gather  ideas  such  as  incorporating  different  service  speeds  to  ship  parts  from  one  echelon 
to  another,  and  the  possible  implementation  of  business  rules  to  determine  which  echelon 
service  customers  require  replacements. 

3.  Simulation 

To  round  out  our  perspective  on  modeling  approaches,  we  examine  a  study  that 
uses  simulation.  Niranjan  and  Ciarallo  formulate  capacitated,  three-echelon,  and  four- 
echelon  systems  with  uncertainty  in  both  demand  and  supply  [17].  Their  approach  uses 
simulation-based  optimization  to  determine  optimal  base  stock  levels  for  the  components 
in  the  system,  based  on  specified  customer  service  levels.  Their  work  gives  us  more 
insight  on  how  to  model  the  evolution  of  inventory  levels  over  time. 
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Overall,  our  work  differs  from  the  cited  literature  in  that  our  formulation 
comprehensively  integrates  uncertainty,  time,  and  monetary  constraints  into  a  single, 
multi-echelon  model,  instead  of  looking  at  these  aspects  in  isolation.  Furthermore,  we 
stray  from  the  literature’s  heavy  reliance  on  the  use  of  expected  values,  in  favor  of 
modeling  more  detailed,  high-fidelity  parameters  and  variables  that  change  over  time. 

C.  PURPOSE 

The  purpose  of  our  research  is  to  develop  an  optimization  model  that  provides 
guidance  on  the  optimal  investment  in  spare  parts  for  a  weapon  system  in  order  to 
maximize  expected  system  availability.  This  includes  modeling  dynamic  interactions  of 
random  part  failures,  time,  and  monetary  constraints  on  the  supply  network.  Even  with 
very  small  networks,  the  number  of  stock  level  decisions  that  must  be  evaluated  ventures 
into  the  tens  of  thousands;  too  many  for  a  human  to  comprehend.  Accordingly,  we  seek 
to  leverage  the  objectivity  and  power  of  mathematical  programming  to  help  decision 
makers  with  this  complex  problem. 

D.  SCOPE,  LIMITATIONS,  AND  ASSUMPTIONS 

1.  Scope 

To  scope  our  model  formulation  and  subsequent  analysis,  we  focus  on  three  key 
areas.  First,  we  only  model  repairable  parts  in  a  closed-supply  network.  This  approach 
aligns  with  a  conclusion  drawn  by  Sherbrooke  that  repairable  parts  most  directly  affect  a 
weapon  system’s  availability,  whereas,  consumable  parts  do  not  [1].  Second,  all  of  our 
echelons  are  in  stationary  locations.  Third,  we  do  not  determine  order  sizes,  reorder 
points,  and  safety  stock  levels  with  classic  supply  management  formulas.  Rather,  we  are 
concerned  with  calculating  optimal  stock  levels  that  result  in  the  best  availability  values 
over  all  randomly-generated  part-failure  scenarios. 

2.  Limitations 

We  limit  the  functionality  of  our  model  to  only  include  the  most  essential  features 
of  a  multi-echelon  supply  network.  As  formulated,  the  model  may  easily  be  modified  to 
include  many  additional  attributes  such  as  more  echelons  and  parts,  implementing  time- 
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varying  transit  times  to  ship  parts  from  one  echelon  to  another,  and  modifying  the 
network  so  that  different  fulfillment  centers  support  different  customers  according  to 
business  rules.  Similarly,  we  have  limited  our  experiments  to  only  a  few  echelons  and 
parts  using  notional  data. 

3.  Assumptions 

To  shape  the  context  of  our  problem  and  accommodate  our  limitations,  the 
following  assumptions  apply  to  our  supply  network’s  policies  and  modeling  logic. 


a.  Policies 

•  Our  budget  only  pays  for  the  initial  number  of  parts  stocked  at  each 
echelon.  This  excludes  inventory  management  expenses  such  as  ordering, 
holding,  and  transportation  costs.  Similarly,  echelons  incur  no  costs  for  the 
future  receipt  of  replacement  parts. 

•  All  failed  parts  are  repaired  by  depots;  none  are  discarded. 

•  Depots  immediately  begin  repair  work  on  failed  parts  that  they  receive 
from  customers. 

•  Repair  and  transit  times  are  deterministic. 

•  Cannibalization  of  weapon  systems  is  not  allowed. 

•  Replacement  parts  for  each  echelon  must  be  satisfied  by  parts  circulating 
within  the  network.  No  external  sources  of  supply  exist. 

b.  Modeling  Logic 

•  If  a  customer  begins  a  time  period  with  insufficient  part(s)  on  hand,  and 
receives  the  necessary  quantity  of  replacements  during  the  same  period, 
the  weapon  system  is  considered  available  during  that  particular  period. 

•  For  the  purposes  of  availability,  part  failures  occur  at  the  end  of  a  time 
period. 

•  A  shipment  supplying  a  particular  replacement  part  to  a  customer  may  be 
initiated  during  the  same  time  period  that  the  installed  part  fails. 

•  A  maximum  of  one  failure  can  occur  per  part,  customer,  and  time  period. 
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II.  METHODOLOGY  AND  FORMULATION 


A.  METHODOLOGY 
1.  Map 

Figure  6  portrays  a  high-level  overview  of  our  research  methodology.  For  any 
weapon  system  we  choose  to  analyze,  we  define  the  identified  sets  and  parameters  as 
inputs,  and  in  return,  our  model  delivers  outputs  that  answer  our  research  questions. 
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Transit  Times 


MODEL 
MAX  Availability 
subject  to: 

Part  Failure  Scenarios 
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Figure  6.  Methodology  map. 


2.  Part  Failures 

Uncertainty  within  the  supply  network  exists  in  many  forms  (e.g.,  part  failure 
uncertainty,  repair  time  uncertainty,  and  transportation  time  uncertainty).  We  focus  our 
efforts  on  modeling  part  failure  uncertainty,  as  part  failures  are  considered  the  most 
uncontrollable  and  detrimental  random  events  that  occur  within  a  supply  network.  When 
a  part  installed  on  a  weapon  system  fails,  the  supply  network  must  respond  to  this  failure 

by  making  a  replacement  part  available  for  issue  as  soon  as  possible,  or  else  availability 
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of  the  weapon  system  suffers.  We  incorporate  part  failures  into  our  formulation  and 
represent  the  evolution  of  our  network  over  time  by  modeling  random  inter-arrival  times 
between  part  failures  for  each  part  installed  at  each  customer.  We  do  this  by  deriving  part 
failure  probability  distributions  from  historical  data  and  using  a  random  part  failure 
generator  to  create  many  randomly  generated  part-failure  scenarios.  Both  are  explained  in 
the  following  section. 

a.  Part  Failure  Probability  Distributions 

Figure  7  displays  example  probability  distributions  of  the  operational  time 
between  failures  for  notional  parts  and  P2  ■  The  horizontal  axes  correspond  to  the 
number  of  time  periods  after  replacement  that  each  part  has  been  operating  and  the 
vertical  axes  give  the  probability  that  each  part  fails  after  operating  for  that  number  of 
time  periods.  Both  parts  have  a  100  percent  probability  of  failing  after  operating  for  35 
time  periods. 


Operational  Time  Periods  Operational  Time  Periods 


Figure  7.  Probability  distributions  of  the  operational  time  between  failures  for  notional 

parts  p\  (left)  and p2  (right). 
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b. 


Random  Part  Failure  Generator 


We  have  developed  a  part  failure  generator  tool  that  accepts  part-failure 
probability  distributions  as  input  and  generates  part-failure  scenarios  as  output.  In 
particular,  the  output  of  this  tool  is  the  parameter/a//_hme^^^;,  which  describes  the 

number  of  time  periods  during  which  echelon  e  has  been  operational  when  part  p  fails  for 
the  time  in  scenario  s.  We  use  this  tool  to  generate  many  scenarios  of  part  failures,  and 
in  turn,  use  the  scenarios  as  input  parameters  for  our  optimization  model.  Our  model  then 
employs  robust  optimization  by  calculating  stock  levels  and  expected  availability  over  all 
input  scenarios.  This  technique  introduces  variation  into  our  metric  calculations  and 
ensures  that  we  do  not  limit  our  results  to  deterministic  solutions  that  are  tied  or 
subordinated  to  a  specific  scenario.  The  overall  goal  of  our  model  is  to  choose  stock 
levels  for  each  echelon  in  such  a  way  as  to  ensure  that  our  supply  chain  performs  well,  on 
average,  over  all  scenarios  considered  simultaneously. 

3.  Part  Failures  and  Availability 

Figure  8  illustrates  the  interactions  between  failures,  transit  times,  and 
availability.  Suppose  we  have  a  weapon  system  comprised  of  parts  and  p^ ,  and  we 
seek  to  calculate  the  system’s  availability  over  the  course  of  eight  time  periods.  We 
assume  that  transit Center, cu.,omer='^^  For  a  particular  scenario  5,  our  random 
part  failure  generator  calculates: 

fail -time 

fail  -time Customer, s=l,i=l=^ 

That  is,  part  p^  is  scheduled  to  fail  for  the  first  time  after  the  system  has  been 

operational  for  two  time  periods,  and  part  p^  is  scheduled  to  fail  for  the  first  time  after 

the  system  has  been  operational  for  three  time  periods.  Accordingly,  we  commence 

operations  and  use  the  system  during  periods  tj  and  .  At  the  end  of  period  >  part  p^ 

fails.  No  replacements  for  part  pi  are  on  hand,  so  the  system  becomes  unavailable  for 

periods  and  while  awaiting  a  replacement  from  the  fulfillment  center.  At  the 
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beginning  of  period  a  replacement  for  pi  arrives  and  is  immediately  installed;  thus, 
the  system  becomes  available  again  during  that  period.  At  the  end  of  period  ,  the 
system  has  been  operational  for  three  time  periods,  so  part  fails.  Again,  no 
replacements  for  part  p^  are  on  hand,  so  transit Center, customer  forces  the  system 
to  become  unavailable  for  time  periods  and  .  At  the  beginning  of  period  ,  a 
replacement  for  pj  arrives  and  is  immediately  installed;  thus,  the  system  becomes 
available  again  during  that  period.  Overall,  the  weapon  system  has  been  available  for  four 
time  periods  out  of  eight.  Per  Equation  1,  the  system’s  availability  is  50  percent. 


Figure  8.  Illustration  of  weapon  system  availability  calculation. 


Note  that  of  our  part  failure  scenarios  do  not  prescribe  the  actual  time  periods  in 
which  failures  occur;  rather,  they  indicate  the  number  of  operational  time  periods  that 
elapse  before  a  part  fails.  The  operational  status  of  the  system  depends,  in  turn,  on  our 
stock  level  decisions.  For  instance,  if  a  spare  of  part  pi  had  been  present  at  time  t2,  our 
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system  would  have  remained  operational  during  period  t^,  and  thus  part  p2  would  have 
failed  at  the  end  of  period  tx  rather  than  at  the  end  of  period  This  subtlety  introduces 
modeling  challenges  not  present  in  most  robust  optimization  problems,  but  it  is  necessary 
in  order  to  accurately  model  the  evolution  of  a  weapon  system. 

B.  FORMULATION 

We  now  describe  our  mathematical  model,  which  uses  stochastic  optimization  to 
maximize  aggregate  weapon  system  availability  by  calculating  optimal  stock  levels  for 
spare  parts  in  a  multi-echelon  supply  network.  It  is  important  to  note  that  our  stock  levels 
account  for  variation  in  the  time  of  failures  because  they  are  chosen  as  the  most  efficient 
values  considering  our  randomly-generated  part  failure  scenarios.  All  other  variables  are 
specific  to  each  scenario. 


1.  Sets  and  Indices 


7?  e  P 
e,e'  eE 
D^E 
F^E 
C^E 
M  ^PxE 

A^ExE 

t^T 
i&I 
s  eS 


part  types 
inventory  echelons 
depots 

fulfillment  centers 
customers 

set  of  part-echelon  pairs  {p,e),  where  echelon  e  can 
process  part  p 

set  of  arcs  (e,e')  joining  echelons  e  and  e'  in  the  shipping 
network 

time  periods 
part  failures 

randomly  generated  scenarios  of  part  failures 


2.  Parameters  [Units] 


repair _timep  ^  number  of  time  periods  required  to  repair  part  p  at  echelon 

e  [time  periods] 

part  _  required minimum  quantity  of  part  p  required  to  operate  a  weapon 
system  at  echelon  e  [parts] 
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transit  _time^y 

costp 

wholesale 

retail 

S 


operational  time  periods  before  part  p  fails  at  echelon  e  for 
the  time  in  scenario  s  [time  periods] 

number  of  time  periods  required  to  ship  a  part  from  echelon 
e  to  echelon  e'  [time  periods] 
cost  of  part  p  [$] 

budget  to  procure  parts  for  depots  and  fulfillment  centers 

[$] 

budget  to  procure  parts  for  customers  [$] 

relative  weight  on  worst-case  scenario  of  availability 
[unitless] 


3.  Decision  Variables  [Units] 


AVAILABLE 

binary;  1  if  weapon  system  is  available  to  operate  at 

echelon  e  during  time  period  t  for  scenario  s,  and  0 
otherwise 

quantity  of  part  p  installed  or  on  hand  at  echelon  e  during 
time  period  t  for  scenario  5  [parts] 

binary;  1  if  INV ^  is  greater  than  or  equal  to 

part  _  required^  ^ ,  and  0  otherwise 

binary;  1  if  fail _time^  ^  ^  ■  equals  'y  AVAILABLE^  ^ ,  and 

t’<t 

0  otherwise 

SHIP,,,.,, 

quantity  of  part  p  shipped  from  echelon  e  to  echelon  e' 
during  time  period  t  for  scenario  5  [parts] 

STOCK  LEVEL^^ 

[j,e 

maximum  quantity  of  part  p  installed  or  on  hand  at  echelon 
e  at  the  beginning  of  the  planning  horizon  [parts] 

WORST  _CASE 

worst-case  system  availability  over  all  scenarios  [fraction] 

4.  Equations 


MAX 

AVAILABLE,  INV , 
ENOUGH,  EAIL,SHIP, 
STOCK _LEVEL,  WORST _CASE 


1  _  o 

,  ,  r,  ,  ,  y  AVAILABLE^,  +5  WORST  CASE  (2) 


subject  to: 
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WORST  CASE  <  — y  AVAILABLE^ , , 

—  'T'\  t,l,S 

^  \  eeC,t 


Vs 


(3) 


INV^,e,^,=  STOCK  _LEVEL^ 


\/(p,e)  ^M,s 


(4) 


e'eC,/l  fail  _timep  g-  ^  j  <t 


repair _time  g— 1,/,5 


e'l(e,e')GA, 

(p,e')£M 


(5) 


\/t  >2,(p,e)  ,s\e  ^  D 


!NV,,,,,,=INV, Z 

e'\{e',e)GA, 

{p,e')GM 


transit  time. 


X  SHIP,,,.,, 


e’\{e,e’)GA, 

{p,e')eM 


(6) 


yt  >2,(p,e)  gM  ,s\e  G  E 


INV^,e,,s  ^  INV^,e,-U  +  Z 

e'\{e'  ,e)GA, 
{p,e')GM 


transit _time^' ^  1  ^^^^p,e,t-\,i,s 

i\fail _timep  ^  ^  ,•  <r-l 

\ft  >  2,{p,e)  eM, 5  I  e  e  C 


(7) 


Z  SHIP,,.,  ,.,  <  X  ™h...o.. 

t'<t, 

i\fail _timep  ^  •  <f' 

V/>,e  e  C,?,^  I  {p,e)  G  M ,  part _  required  ^  >  1 


t'<t,e'\{e',e)GA 

{p,e')eM 


p,e 


(8) 


Z  SHIP,,,,,<  X  FAIL,,.,^ 


t'<t, 

e'l(e',e)eA, 

(p,e')eM 


t'<t, 

e'eCl(e,e')£A, 
i\fail _time  p  ^  j  <t' 


Z  ship,,,,.,<  Y.  fail,,,.,. 


t'<t, 

e'\{e,e')GA, 

(p,e')GM 


t'<t,e'GC, 
ilfail  _time  p  ^  j  <t' 


yp,e  GE,t,s\  ip,e)  gM 


yp,e  ^D,t,s\  (p,e)  gM 


(9) 


(10) 


AVAILABLE^  ,^^<  ENOUGH 

\/p,  e^C,t,s\  (p,  e)  G  M ,  part  _  required  p  ^  >  1 


(11) 


AVAILABLE > 


ENOUGH 


p\{p,e)GM , part _  required  ^  >1 


p,e,t,s 


Z  1 

\^p\{p,e)GM ,part _requiredp  g>l  J 

\/e  GC,t,s 


+  1 


(12) 
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INV 

ENOUGH 


p,e,t,s 


part _required p  ^ 

V77,  e  ^C,t,s\  {p,e)  G  M , part _required p  ^  >  1 


(13) 


ENOUGH > 


INVp  e  r  s  “  part  _  required ^  ^  +  0.99 


1  +  part  _  required p  ^  + 


retail 
cost  „ 


'sfp,e  ^C,t,s\{p,e)  G  M ,  part  _  required  p  ^  >  1 


(14) 


fail  _ time p “  '^AVAILABLE ^ ,, 


t'<t 


\/p,e  G  C,t,i,s\  ip,e)  & M , part _ required ^  ^  >  I,  fail _time ^  e  .5,  -  ^ 


p,e,s,i 


AVAILABLE,,, 


fAIL,„j,  < 


t'<t 


faU_time^^^. 

\/p,  e  GC,t,i,s\(p,e)  gM  ,  part  _  required  ^ fad  _  time^  ^  ^ ,  <t 


Z  fail,,,,,> 

t'\fail  _timep  ^  ^  j  <t'<t 


1  -  fail _  time^  ,,  ,,j  +  ^AVAILABLE,,  , 


t'<t 


t  -  fail  _  time  ^ + 1 


yp,e,i,s  \{p,e)  ^ M , part _required p ,,  >  1 


(15) 


(16) 


^  FA/L^^,.^<1  yp,e,i,s\{p,e)GM,part_requiredp,^>l  (17) 

t>fail _time  ^  ^  j 


(18) 


INV STOCK _LEVEL^^ 


\/p,e,t,s\{p,e)  &M  ,e  (19) 


STOCK  _  LEVEL ^  ^  >  part  _  required ^  ^ 

\/p,e  gC,s  \{p,e)  G  M , part _required p  ^  >  1 

V  STOCK  LEVEL,,  Nost,  <  wholesale 

{p,e)^M\eiC 


(20) 

(21) 


^  STOCK _LEVELp,,costp<  retail  +  ^  part _required ^  ,,cost ^  (22) 


(p,e)EMleGC, 


(p,e)eM  leeC, 
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,  ,  e{0,l} 

\/e,t,s 

(23) 

ENOUGH e{0,l} 

Vp,e,t,s 

(24) 

yp,e,t,s 

(25) 

\/p,e,t,i,s 

(26) 

\/p,e,t,s 

(27) 

STOCK  _  LEVEL^^  >0  \/p,e 

(28) 

WORST  _CASE>0 

(29) 

5.  Objective  Function 

Our  objective  in  Equation  2  is  to  maximize  the  weighted  average  of  the 
availability  of  our  weapon  system  over  all  part-failure  scenarios  and  availability  for  the 
worst-case  scenario.  As  calculated  in  Equation  3,  use  of  WORST _CASE  allows  us  to 
limit  overestimation  of  availability  over  extreme  values.  In  our  computational  examples, 
we  typically  set  5  equal  to  0. 1 . 

6.  Constraints 

a.  Inventory  Levels 

Equation  4  initializes  inventory  levels  of  each  echelon  in  the  first  time  period.  In 
accordance  with  Eigure  8,  Equations  5  through  7  use  fail  _time  ^  ^  ^  ^ ,  transit  _time^  ,  and 

repair _time p  ^ ,  to  model  the  flow  of  parts  through  our  supply  network.  Equation  5 

calculates  current  inventory  levels  for  depots  as  failed  parts  are  received  from  customers, 
repaired,  and  issued  to  fulfillment  centers  as  replacements.  Similarly,  Equation  6 
calculates  current  inventory  levels  for  fulfillment  centers  as  parts  are  received  from 
depots  and  subsequently  issued  as  replacements  to  customers.  Einally,  Equation  7 
calculates  current  inventory  levels  for  customers  as  part  failures  occur  and  replacements 
are  received  from  supporting  fulfillment  centers. 
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b.  Shipping 

Equations  8  through  10  prevent  our  model  from  anticipating  and  preemptively 
responding  to  part  failures  that  occur  in  future  time  periods.  This  guarantees  that 
prescient  decisions  do  not  falsely  inflate  our  availability  values  based  on  knowledge  we 
would  not  have  in  a  real  situation.  Specifically,  Equation  8  ensures  that  no  customer 
receives  more  of  part  p  than  the  total  number  of  failures  they  have  had  for  part  p  as  of 
time  t.  Equation  9  specifies  that  no  fulfillment  center  can  receive  more  of  part  p  than  the 
sum  of  failures  their  customers  have  had  for  p  as  of  time  t.  Einally,  Equation  10  states 
that  the  total  quantity  of  shipments  of  part  p  out  of  all  depots  cannot  exceed  the  total 
number  of  failures  of  p  over  all  customers  as  of  time  t. 

c.  Weapon  System  Availability 

Equations  11  through  14  determine  whether  our  weapon  systems  have  sufficient 
quantities  of  each  part  on  hand  for  each  customer.  Specifically,  in  Equation  13,  if  the 
current  inventory  of  part  p  is  less  than  the  quantity  of  p  required  to  operate  the  weapon 
system,  the  ENOUGH  variable  records  this  shortfall.  Equations  11  and  12  ensure  that  the 
system  becomes  unavailable  until  enough  parts  are  on  hand. 

d.  Failure  Events 

Equations  15  through  18  signal  the  time  periods  when  failures  occur  for  each  part 
in  operation  at  each  customer  echelon.  The  complexity  of  these  equations  involves 
ensuring  that  each  part  fails  in  accordance  with  our  randomly-generated  part  failure 
scenarios.  We  accomplish  this  with  the  parameter  fail  _time^  ^  ^  ■  and  the  binary  variables 

AVAILABLE^  ^  ^  and  EAIL^  ^  ^ ^ .  Equation  15  prevents  part  failures  from  occurring  when 
the  sum  of  operational  time  periods  up  to  time  t  is  strictly  greater  than  fail  _time ^  ^  ^  . 
Conversely,  Equation  16  prevents  part  failures  from  occurring  when  the  sum  of 
operational  time  periods  up  to  time  t  is  strictly  less  than  fail_timep^^..  Equation  17 
ensures  that  part  p  at  echelon  e  cannot  fail  more  than  once  for  the  time  over  all  time 
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periods  in  each  scenario.  Finally,  Equation  18  forces  ^  to  signal  a  part  failure 

when  fail  _time^  ^  ^  ■  equals  the  sum  of  operational  time  periods  up  to  time  t. 

e.  Stock  Levels 

Equations  19  and  20  limit  the  values  of  our  STOCK _LEVEL^  ^  decision  variable. 
Specifically,  Equation  19  ensures  that  never  exceeds  STOCK  _LEVEL^^, 

regardless  of  the  scenario.  Eurthermore,  part  _required  ^  ^specifies  the  quantity  of  each 

part  that  must  be  on  hand  in  order  for  a  customer  to  operate  its  weapon  system.  Similar  to 
Equation  19,  we  use  Equation  20  to  ensure  that  STOCK _LEVELp ^  is  always  greater 
than  or  equal  to  this  quantity. 

f  Budget 

Equations  21  and  22  set  wholesale  and  retail  budget  limits. 
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III.  ANALYSIS 


To  conduct  our  analysis,  we  implement  our  formulation  using  the  Generalized 
Algebraic  Modeling  System  (GAMS)  version  23.8.2  [18]  and  solve  it  as  a  mixed-integer 
linear  program  with  the  CPLEX  12.4.00  solver  [19]. 

A.  RESULTS 

1.  Investment  versus  Availability 

To  gain  insights  similar  to  those  discussed  in  Section  I.A.2.C,  we  temporarily 
simplify  our  formulation  and  consider  a  single  budget  constraint  instead  of  separate 
wholesale  and  retail  budgets.  That  is,  we  combine  Equations  21  and  22  and  limit  the  total 
cost  to  purchase  spare  parts  for  all  echelons  by  a  given  budget  value.  Table  1  displays 
part  attributes  we  use  as  model  inputs  for  our  experiments,  and  we  model  the  flow  of 
these  parts  throughout  the  supply  network  shown  in  Eigure  10(d),  over  15  time  periods 
and  15  scenarios.  Each  customer  requires  one  of  each  part  to  operate  their  weapon 
system,  and  our  investment  budget  ranges  from  $0  to  $100  in  $5  increments.  Eigure  9 
illustrates  the  resulting  investment  versus  availability  curve. 


Part  Types 

Average  Failure  Time 
Poisson  Distribution 
(Time  Periods) 

repair _timep  Depot 
(Time  Periods) 

COStp 

($) 

Pi 

1 

1 

6 

P2 

2 

2 

1 

P3 

3 

3 

3 

Table  1.  Repairable  part  attributes  for  experiments. 
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Figure  9.  Investment  versus  availability  curve. 


Intuitively,  our  model  produces  results  as  expected.  As  our  investment  budget 
increases,  our  model  chooses  to  purchase  more  spare  parts  to  stock  at  echelons  that  result 
in  the  best  availability  values.  Accordingly,  we  have  more  spares  on  hand  to 
accommodate  part  failures,  so  availability  ultimately  improves.  Clearly,  our  model 
exhibits  the  investment  versus  availability  relationship  discussed  in  Chapter  I.  Note  that 
if  we  choose  to  invest  $0  into  purchasing  spare  parts,  we  estimate  that  our  system  will  be 
available  33  percent  of  the  time.  Conversely,  if  we  choose  to  invest  $85  optimally,  we 
estimate  that  our  system  will  be  available  100  percent  of  the  time.  Any  investment 
beyond  $85  does  not  increase  availability  for  the  scenarios  considered. 

Note  that  increasing  the  budget  does  not  always  result  in  availability  gains.  Recall 
that  our  formulation  is  based  on  the  integrated  modeling  of  uncertainty,  time,  and 
monetary  constraints  on  a  multi-echelon  supply  network’s  ability  to  support  a  weapon 
system.  During  any  given  time  period,  failures  occur,  parts  are  repaired,  replacements  are 
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in  transit,  and  our  model  must  choose  the  best  stock  levels  that  accommodate  all  of  these 
events  over  all  time  periods  and  scenarios.  Purchasing  many  parts  of  one  type  may  not 
result  in  as  great  of  an  impact  on  availability  as  purchasing  one  of  another  type  of  part. 
For  each  $5  increment  in  our  investment  budget,  part  types  compete  against  each  other. 
We  purposely  stagger  our  part  costs  from  $1  to  $6  to  exaggerate  this  effect  and  ensure 
that  our  model  cannot  buy  one  more  of  each  part  with  each  budget  increment.  Thus, 
availability  gains  are  not  guaranteed  each  time  the  investment  budget  increases. 

2.  Multi-Echelon  Supply  Network  Experiments 

Next,  we  experiment  with  different  sizes  of  multi-echelon  supply  networks  to 
realize  the  impacts  of  topology  and  part  types  on  availability.  To  do  this,  we 
reincorporate  our  wholesale  and  retail  budgets  and  model  part  failure  scenarios  over  each 
network  shown  in  Figure  10.  Table  2  itemizes  initial  conditions  for  our  experiments.  In 
order  to  reduce  computational  time,  we  increase  the  absolute  tolerance  {availability  gap) 
for  larger  networks.  Each  customer  is  assumed  to  require  one  of  each  part  to  operate  their 
weapon  system. 


29 


o, 
\  * 

II  s 


Dqwt 

1 


Fulfillmenl  Center 


Castomcr: 


(a) 


transit  time=  1 


Depot 


Fullillmeoi  Center  1 


Custotner  1 


FuUilhtkent  Center  2 


(b) 


V 


CusionKr2 


transit  time=  I 


(c) 


transit  time-  1 


(d) 


transit  time  -  1 


—  —  —  transit  time  =  2 


—  —  —  transit  time  =  2 


Figure  10.  Multi-echelon  supply  networks  for  experiments. 


Network 

Topology 

s 

Time  Periods 

Scenarios 

Absolute 

Tolerance 

(%) 

Four-Echelon 

Figure  10(a) 

0.1 

15 

15 

5 

Five-Echelon 

Figure  10(b) 

0.1 

15 

15 

5 

Seven-Echelon 

Figure  10(c) 

0.1 

15 

15 

10 

Eight-Echelon 

Figure  10(d) 

0.1 

15 

15 

15 

Table  2.  Initial  conditions  for  experiments. 
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Figure  11  displays  our  estimated  availability  values  for  each  combination  of 
wholesale  and  retail  budgets  for  weapon  systems  comprised  of  parts  pi  and  p2.  Likewise, 
Figure  12  displays  availability  for  weapon  systems  comprised  of  parts  pu  p2,  and  ps.  The 
horizontal  axes  indicate  our  retail  budget  and  vertical  axes  indicate  our  wholesale  budget. 
Availability  values  are  shown  at  the  intersection  of  each  wholesale  and  retail  budget 
combination. 

Looking  at  Figure  11,  we  conclude  that  as  supply  network  topology  becomes 
more  complex,  more  investments  in  spare  parts  are  needed  to  ensure  high  rates  of 
availability.  Furthermore,  we  see  that  regardless  of  network  size,  if  we  choose  to  invest 
$0  into  purchasing  spares  for  any  echelon,  our  weapon  system  availability  will  never 
exceed  42  percent.  The  plots  give  us  insight  into  the  importance  of  having  a  sufficient 
retail  budget.  For  example,  if  we  choose  to  invest  $0  into  the  purchase  of  spares  for  all 
customers,  no  matter  how  much  we  invest  into  stocking  spares  at  depots  and  fulfillment 
centers,  we  estimate  availability  will  never  exceed  79  percent.  This  is  because  weapon 
systems  requiring  replacement  parts  will  always  be  unavailable  for  at  least 

transit  _  time center, customer  ' 
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Color  Key 

EJ 

0.4  0.6  0.8  1 
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Figure  11. 


Availability  estimates  modeling  parts  pi  and  p2. 
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Figure  12.  Availability  estimates  modeling  parts  pi,  p2,  and  ps. 


Our  experiments  with  three  part  types  also  show  that,  as  the  network  beeomes 
more  complex,  more  investment  is  needed  to  achieve  a  certain  availability.  Moreover, 
regardless  of  network  size,  if  we  choose  to  invest  $0  into  purchasing  spares  for  any 
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echelon,  at  best,  our  weapon  system  availability  will  never  exceed  38  percent.  Comparing 
Figures  1 1  and  12,  we  see  how  modeling  more  part  types  results  in  visible  trends  of  lower 
availability. 

The  general  insight  displayed  by  all  of  these  networks  is  that  for  any  fixed 
wholesale  budget,  availability  increases  at  a  faster  rate  when  we  increase  our  retail 
budget  than  when  we  increase  our  wholesale  budget.  We  reason  that  this  is  most  heavily 
influenced  by  transit  cemer. customer-  Recall  that  when  customers  have  spare  parts 

on  hand,  they  can  immediately  replace  failed  parts  and  keep  their  weapon  system 
operational,  as  opposed  to  waiting  for  replacement  parts  to  arrive  from  a  fulfillment 
center.  Even  though  a  customer’s  availability  suffers  while  the  part  is  in  transit,  a  proper 
balance  of  spare  part  stock  must  be  struck  between  depots  and  fulfillment  centers  because 
they  have  the  ability  to  supply  multiple  customers. 

To  highlight  efficiencies  gained  from  multi-echelon  modeling,  we  explore 
availability  values  from  Figure  12(d),  which  corresponds  to  an  eight-echelon  network 
modeling  the  flow  of  parts  pi,  p2,  and  ps.  Suppose  we  wish  to  partition  a  given  total 
budget  optimally  into  wholesale  and  retail  budgets.  From  Figure  12(d),  we  note  that  each 
partition  may  result  in  a  different  availability  value.  Figure  13  illustrates  an  investment 
versus  availability  curve  where  each  data  point  represents  the  availability  achieved  for  a 
particular  partition  of  the  budget  shown  on  the  horizontal  axis. 
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Figure  13.  Investment  versus  availability  eurve  for  an  eight-echelon  network  modeling 
parts  pi,  p2,  and  ps-  Each  data  point  represents  a  particular  division  of  the 
total  budget  indicated  on  the  horizontal  axis  into  wholesale  and  retail 

components. 

Table  3  lists  the  particular  availability  values  that  we  examine.  From  the  data,  we 
see  that  if  a  decision  maker  chooses  to  optimally  invest  any  budget  between  $80  and  $95 
into  purchasing  spares,  weapon  system  availability  is  expected  to  be  96  percent.  Hence,  a 
prudent  choice  is  to  only  invest  $80  and  use  the  $15  savings  elsewhere.  Furthermore,  if 
the  decision  maker  chooses  to  accept  a  one  percent  reduction  in  availability  from  96 
percent  to  95  percent,  an  additional  savings  of  $5  can  be  realized.  Investing  $70  or  less, 
results  in  a  large  drop  in  availability.  We  generalize  that  by  taking  advantage  of 
efficiencies  gained  from  multi-echelon  modeling,  for  a  smaller  investment  in  spares,  we 
may  obtain  the  same  or  a  comparable  level  of  weapon  system  effectiveness. 
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$15  savings  for 
same  Availability 


Total 

Investment 

($) 

WTiolesale 

Budget 

($) 

Retail 

Budget 

($) 

Availability 

(%) 

95 

45 

50 

96 

90 

45 

45 

96 

85 

35 

50 

96 

80 

35 

45 

96 

75 

25 

50 

95 

70 

40 

30 

89 

$20  savings  if  the  decision 
maker  accepts  a  1% 
reduction  in  Availability 


Table  3.  Availability  for  an  eight-echelon  network  modeling  parts  pu  p2,  and  py. 


B.  SOLVE  TIME 

Because  we  choose  to  model  many  of  the  detailed  interactions  that  take  place 
within  a  multi-echelon  supply  network,  our  formulation  is  challenged  by  lengthy  solve 
time.  To  gain  an  appreciation  for  this,  Table  4  displays  the  problem  size  of  our 
experiments.  Solve  times  for  these  instances  are  on  the  order  of  days.  Clearly,  even  for 
our  small  experiments,  finding  an  optimal  solution  is  a  substantial  undertaking  for  the 
solver. 


Network 

Topology 

s 

Part 

Types 

Variables 

Scenarios 

Discrete 

Variables 

Constraints 

Four- 

Echelon 

Figure  10(a) 

0.1 

Pb  P2 

7,399 

15 

5,328 

15,167 

0.1 

Pi,  P2,  P3 

10,205 

15 

7,273 

20,723 

Five- 

Echelon 

Figure  10(b) 

0.1 

Pi,  P2 

9,141 

15 

6,678 

18,317 

0.1 

Pi,  P2,  P3 

12,908 

15 

9,298 

25,448 

Seven- 

Echelon 

Figure  10(c) 

0.1 

Pi,  P2 

15,180 

15 

11,626 

33,196 

0.1 

Pi,  P2,  P3 

20,990 

15 

14,394 

45,207 

Eight- 

Echelon 

Figure  10(d) 

0.1 

Pi,  P2 

15,180 

15 

10,992 

33,196 

0.1 

Pi,  P2,  P3 

20,990 

15 

14,394 

45,207 

Table  4.  Problem  size  of  multi-echelon  experiments. 


To  overcome  this  challenge,  we  examine  our  formulation  for  improvements  that 
may  reduce  solve  time.  The  user  specifies  most  of  the  input  data  (e.g.,  part  types,  part 
failure  probabilities,  network  topology,  time  periods,  and  budgets),  so  we  do  not  consider 
modifying  those  items.  However,  an  aspect  that  we  can  easily  change  is  the  number  of 
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part-failure  scenarios  over  which  we  optimize.  Conceptually,  the  more  scenarios  we 
consider  simultaneously,  the  more  accurate  our  availability  estimate  will  be;  however, 
more  scenarios  result  in  longer  solve  times.  Accordingly,  we  seek  to  determine  an 
adequate  number  of  scenarios  that  strikes  a  reasonable  balance  between  availability 
precision  and  solve  time. 

C.  AVAILABILITY  AND  SOLVE  TIME  VARIATION 

To  gain  an  understanding  of  the  variation  in  our  output,  we  plot  confidence 
intervals  for  availability  and  solve  time  with  respect  to  the  number  of  part  failure 
scenarios  over  which  we  optimize.  We  accomplish  this  by  conducting  independent  (i.e., 
different)  replications  for  each  experiment. 

We  perform  this  analysis  with  the  four-echelon  supply  network  shown  in  Figure 
10(a)  and  model  the  flow  of  parts  pi  and  p2  over  a  horizon  of  15  time  periods.  Table  5 
displays  part  attributes.  Wholesale  and  retail  budgets  are  fixed  at  $5  each,  and  we  solve  to 
an  absolute  tolerance  of  five  percent  availability.  Furthermore,  we  run  100  replications 
for  each  scenario,  and  our  scenario  counts  range  from  1  to  50  in  increments  of  five. 
Figure  14  plots  95  percent  confidence  intervals  for  availability  and  solve  time,  given  a 
particular  number  of  scenarios  for  each  replication.  The  horizontal  axes  indicate  the 
number  of  part  failure  scenarios  used  within  each  replication  and  the  vertical  axes 
indicate  availability  and  solve  time,  respectively. 


Average  Failure  Time 

repair_timep  Depot 

Part  Types 

Poisson  Distribution 

COStp 

(Time  Periods) 

(Time  Periods) 

($) 

Pi 

1 

6 

6 

P2 

2 

3 

1 

Table  5.  Repairable  part  attributes  for  variation  experiment. 
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Scenano  Count 


(a) 


Scenano  Count 

(b) 

Figure  14.  95  percent  confidence  intervals  for  availability  and  solve  time. 

With  respect  to  availability,  we  note  that  as  the  number  of  part  failure  scenarios 
used  increases  from  1  to  50,  our  confidence  intervals  narrow  and  our  availability 
estimates  become  more  precise.  Noticeably,  if  we  choose  to  model  only  one  scenario,  we 
expect  our  results  will  have  a  large  range  of  variation  and  our  availability  estimates  will 
be  overly  optimistic  (i.e.,  subordinating  our  stock  level  decisions  to  that  solution  will  be 
highly  suboptimal).  The  confidence  interval  for  one  scenario  is  over  twice  the  width  of 
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that  of  five  scenarios,  yet  the  computational  time  to  solve  is  approximately  the  same. 
Thus,  optimizing  over  multiple  scenarios  simultaneously  is  recommended.  Our  mean 
availability  estimates  and  confidence  interval  widths  reach  a  steady  state  when  we 
employ  25  or  more  scenarios  per  replication.  This  informs  us  that  modeling  more  than  25 
scenarios  would  be  an  unnecessary  use  of  computing  resources. 

In  stark  contrast  to  our  availability  confidence  interval  plots,  we  see  that  the 
trajectory  of  the  average  solve  time  per  replication  resembles  that  of  an  exponential 
curve.  Furthermore,  as  the  scenario  count  per  replication  increases,  so  does  the  variation 
in  the  amount  of  time  it  takes  to  solve  each  replication.  A  drastic  change  in  solve  time 
variation  occurs  between  25  and  30  replications:  the  confidence  interval  for  30  scenarios 
per  replication  is  nearly  four  times  as  wide  as  the  interval  for  25  scenarios  per  replication. 
Combining  this  observation  with  the  fact  that  our  availability  confidence  intervals  reach  a 
steady  state  at  25  scenarios  per  replication,  we  reason  that  modeling  25  scenarios  per 
replication  gives  us  outputs  that  reasonably  balances  availability  precision  and  solve  time 
for  this  particular  experiment. 

We  generalize  these  two  notions  and  conclude  that  for  any  particular  supply 
network,  part  types,  time  horizon,  and  budget  combination  we  choose  to  model,  there 
exists  a  scenario  count  that  strikes  a  reasonable  empirical  balance  between  availability 
precision  and  solve  time,  making  the  modeling  of  more  scenarios  unnecessary. 

D.  STOCK  LEVEL  CONVERGENCE 

In  addition  to  balancing  availability  precision  and  solve  time,  we  predict  that  the 
number  of  part  failure  scenarios  used  also  plays  a  role  in  our  convergence  to  an  optimal 
solution  for  STOCK _LEVEL^ Conceptually,  we  reason  that  as  scenario  count 

increases,  our  model  should  more  frequently  choose  a  specific  stock  level  plan  for  each 
spare  part  at  each  echelon.  To  illustrate  this  concept,  we  conduct  an  experiment  with  the 
four-echelon  supply  network  shown  in  Figure  10(a)  and  model  the  flow  of  parts  pi  and  p2 
with  wholesale  and  retail  budgets  fixed  at  $5  each.  We  run  100  replications  for  each 
scenario  where  our  scenario  counts  range  from  1  to  50  in  increments  of  five,  and  we 
solve  to  an  absolute  tolerance  of  one  percent  availability.  Figures  15  and  16  display  the 
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resulting  wholesale  and  retail  stock  level  histograms  and  corresponding  plans.  The 
horizontal  axes  indicate  the  number  of  part  failure  scenarios  used  within  each  replication 
and  the  vertical  axes  indicate  the  number  of  times  a  particular  stock  level  plan  is  chosen 
as  the  optimal  solution  for  each  replication.  Stock  level  plans  accompany  each  histogram. 
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Figure  15.  Wholesale  stock  level  histogram  and  plans. 
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Figure  16.  Retail  stock  level  histogram  and  plans. 


Per  the  figures,  we  can  see  that  as  scenario  count  increases,  our  model  more 
frequently  chooses  a  particular  stock  level  plan  (i.e.,  converges  towards  an  optimal 
solution).  Plans  W-A  and  R-A  are  most  frequently  chosen  as  the  plans  that  result  in  the 
best  availability  values  when  we  optimize  over  one  part  failure  scenario.  However,  recall 
that  we  anticipate  it  may  be  highly  suboptimal  to  optimize  over  one  scenario  only;  so,  we 
focus  our  attention  on  results  where  scenario  count  is  greater  than  one.  Note  that  as  we 
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increase  the  scenario  count  from  10  to  50,  plans  W-A  and  R-A  steadily  emerge  as  the 
most  frequently  chosen  solutions  while  other  stock  level  plans  subside.  Thus,  we 
conclude  that  for  this  particular  experiment,  our  model  converges  on  stock  level  plans  W- 
A  and  R-A  as  the  optimal  solutions  for  STOCK _LEVEL^ Note  that  for  the 

experimental  setup  we  have  considered  multiple  optimal  solutions  are  possible. 

E.  VALUE  OF  STOCHASTIC  MODELING 

To  demonstrate  the  value  of  stochastic  modeling,  we  now  consider  a  simple 
system  in  which  only  10  part  failure  scenarios  are  possible.  In  this  simplified  setting,  we 
compare  availability  values  resulting  from  optimization  over  a  single  scenario  vice 
optimization  over  all  10  scenarios.  We  model  the  flow  of  parts  pi,  p2,  and  ps  throughout 
the  supply  network  shown  in  Figure  10(b)  over  15  time  periods  and  solve  to  an  absolute 
tolerance  of  five  percent  availability.  Table  6  displays  part  attributes  for  our  comparisons, 
and  wholesale  and  retail  budgets  are  fixed  at  $10  and  $15,  respectively. 


Part  Types 

Average  Failure  Time 
Poisson  Distribution 
(Time  Periods) 

repair _timep  Depot 
(Time  Periods) 

COStp 

($) 

Pi 

1 

6 

6 

P2 

2 

3 

1 

P3 

3 

1 

3 

Table  6.  Repairable  part  attributes  for  stochastic  modeling  experiment. 


Let  5  and  t  denote  part  failure  scenarios  from  the  set  5  =  |5'j,5'2,...,5'jo} . 

Furthermore,  let  STOCK  _  LEVEK^  ^  denote  the  optimal  stock  level  plan  resulting  from 

optimizing  only  over  scenario  5,  and  let  STOCK _LEVEL^^  denote  the  optimal  stock 

level  plan  resulting  from  optimizing  over  all  10  scenarios  in  S.  If  we  optimize  over  only  a 
single  scenario  s  and  implement  the  resulting  stock  level  plan,  we  cannot  be  certain  that  s 
is  sufficiently  representative  of  all  possible  scenarios.  Should  a  different  scenario  occur, 
we  may  experience  a  substantially  different  availability  value  than  anticipated,  and  one 
that  is  much  lower  than  could  have  been  achieved  with  perfect  knowledge.  To  quantify 
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this  difference,  let  availability ^  ,  denote  the  availability  value  that  results  when  we 
implement  stock  level  plan  STOCK _LEVEIl^  ^ ,  but  scenario  5  actually  occurs.  Likewise, 
let  availability^  denote  the  availability  achieved  when  we  implement 
STOCK  _LEVEL^^^  and  scenario  5  actually  occurs. 

Figure  17  displays  the  results  of  our  ten- scenario  experiment.  The  horizontal  axis 
indicates  the  scenario  5  that  actually  occurred,  while  the  data  points  indicate 
availability^ ,  for  all  t  e  5 ,  as  well  as  availability^  .  As  the  figure  indicates, 

availability ^  is  relatively  high  for  all  scenarios;  in  fact,  for  many  scenarios, 

availability^  =  availability^  ^ .  Furthermore,  the  risk  of  achieving  a  highly  suboptimal 
availability  value  is  much  lower  when  using  STOCK  _LEVEL^^^^ .  As  the  figure 

indicates,  for  a  single  outcome  5,  we  see  an  absolute  difference  in  availability  of  up  to  30 
percent  depending  on  which  scenario  t  was  used  in  the  optimization;  in  contrast,  we  have 
at  most  a  14  percent  absolute  difference  between  availability ^  and  availability^  ^  for 
all  scenarios.  Table  7  summarizes  the  results  of  our  experiment. 


43 


o 


Scenario 


Figure  17.  Availability  plot  for  stochastic  modeling  experiment. 


Scenario 

s 

availability^  ^ 

Worst  Case 
availability^ , 

Average 
availability ^ , 

availability^ 

Si 

0.97 

0.67 

0.75 

0.83 

S2 

0.93 

0.67 

0.79 

0.93 

S3 

0.83 

0.63 

0.73 

0.83 

S4 

0.77 

0.60 

0.71 

0.77 

S5 

0.77 

0.67 

0.72 

0.77 

S6 

0.73 

0.57 

0.66 

0.73 

s? 

0.70 

0.63 

0.67 

0.70 

Ss 

0.70 

0.53 

0.61 

0.60 

Sg 

0.63 

0.53 

0.58 

0.60 

SlO 

0.63 

0.47 

0.55 

0.53 

Table  7.  Availability  results  for  stochastic  modeling  experiment. 
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F.  TECHNIQUES  TO  REDUCE  SOUVE  TIME 

Because  our  formulation  models  detailed  interactions  within  a  multi-echelon 
supply  network,  we  face  challenges  in  terms  of  solve  time.  We  implement  restriction  and 
relaxation  techniques  to  help  our  solver  find  an  optimal  solution  faster. 

For  a  restriction,  we  temporarily  fix  our  STOCK  _LEVEL^  ^  decision  variables  to 

obtain  an  initial  solution,  and  then  use  the  corresponding  variable  bounds  as  an  input,  to 
obtain  a  full  solution  for  our  mixed- integer  program.  We  accomplish  this  by  fixing 
STOCK _LEVELp  ^  to  equal  the  maximum  quantity  of  spares  we  could  purchase  for  each 

part  at  each  echelon,  given  our  wholesale  and  retail  budgets.  By  solving  this  restriction, 
we  obtain  an  initial  integer  solution  for  our  problem,  which  occasionally  helps  the 
solver’s  optimization  process. 

In  terms  of  implementing  a  relaxation,  we  suspend  our  shipping  constraints  to 
establish  an  upper  bound  on  our  optimal  objective  value.  Then,  we  use  this  upper  bound 
to  incrementally  refine  bounds  on  our  decision  variables  so  that  the  solver  more  quickly 
converges  on  an  optimal  solution.  We  accomplish  this  by  conducting  a  preliminary  solve 
that  excludes  Equations  8  through  10.  We  reason  that  by  suspending  the  shipping 
constraints,  we  create  an  environment  where  we  foresee  exactly  when  and  where  part 
failures  will  occur.  In  response,  our  model  takes  preemptive  action  and  ships  replacement 
parts  to  customers  where  failures  are  about  to  occur.  Conceptually,  our  availability 
estimate  should  never  be  greater  than  in  situations  where  we  can  see  the  future  and  take 
action  accordingly.  The  discussed  approaches  result  in  modest  reductions  to  solve  time. 
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IV.  SUMMARY  AND  FUTURE  RESEARCH 


A.  SUMMARY 

When  many  identical  weapon  systems  operate  in  parallel  and  rely  on  a  multi¬ 
echelon  supply  network  for  replacement  parts,  decision  makers  face  the  problem  of 
determining  how  to  optimally  choose  stock  levels,  given  a  limited  budget. 

Our  research  confronts  this  problem  by  formulating  a  mathematical  model  that 
determines  the  optimal  investment  of  limited  budgets  into  the  procurement  of  spare  parts 
so  that  weapon  system  availability  is  maximized.  Specifically,  we  use  stochastic 
optimization  to  integrate  the  modeling  of  uncertainty,  time,  and  monetary  constraints 
onto  a  multi-echelon  supply  network.  First,  we  incorporate  uncertainty  into  our  model  by 
randomly  generating  scenarios  of  time-based  part  failures  (for  a  particular  weapon 
system),  which  become  input  parameters  for  our  stock  level  calculations.  Second,  we 
model  the  evolution  of  our  inventory  levels  over  time  and  incur  a  penalty  in  our  objective 
function  if  the  system  is  unavailable  (i.e.,  if  the  supply  network  does  not  provide 
replacement  parts  in  a  timely  manner).  Accordingly,  we  provide  decision  makers  with  an 
estimate  of  expected  availability  given  an  optimal  investment  in  spare  parts  for  the 
system. 

Rather  than  looking  at  supply  network  features  such  as  part  failures,  repair,  transit 
times,  and  costs  in  isolation,  our  formulation  integrates  their  interactions  in  an 
optimization  model  that  seeks  maximizing  overall  expected  availability.  Our 
computational  experiments  show  that,  as  supply  network  topology  becomes  more 
complex,  more  investments  in  spare  parts  are  needed  to  ensure  high  rates  of  weapon 
system  availability.  We  show  the  tradeoff  relationship  between  solution  accuracy  and  the 
time  it  takes  to  find  an  optimal  solution  as  the  number  of  scenarios  considered  increases. 

B.  FUTURE  RESEARCH 

We  view  our  research  as  the  first  step  towards  detailed  multi-echelon  modeling 
that  is  not  otherwise  examined  in  the  literature.  Future  research  may  build  upon  our  work 

in  terms  of  algorithmic  refinements,  technical  enhancements,  and  further  analysis. 
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1.  Algorithmic  Refinements 

We  propose  that  future  work  build  upon  our  efforts  to  decrease  solve  time. 
Although  we  have  performed  a  preliminary  investigation  of  techniques  to  improve  solve 
time,  we  believe  that  more  can  be  done  to  improve  the  model’s  efficiency.  Possible 
refinements  may  include  reformulation  of  some  constraints,  especially  those  relating 
failures  and  availability,  further  experimentation  with  heuristic  and  exact  solution 
methods,  as  well  as  decomposition  approaches. 

2.  Technical  Enhancements 

Given  that  we  seek  to  accurately  estimate  availability,  we  envision  the  addition  of 
time-varying  transit  times  and  mobile  customers  as  future  technical  enhancements  to  our 
model.  As  it  stands,  our  model  only  accounts  for  part  failure  uncertainty  and  neglects  any 
randomness  associated  with  the  number  of  time  periods  it  takes  to  ship  a  replacement  part 
from  one  echelon  to  another.  In  order  to  better  estimate  weapon  system  availability,  we 
must  account  for  this  uncertainty.  We  foresee  that  the  scope  of  our  part  failure  scenarios 
should  expand  to  include  randomly  generated  transit  times. 

Additionally,  in  our  current  model,  each  customer  is  supported  by  the  same 
fulfillment  center  over  all  time  periods.  This  assumption  oversimplifies  reality  because 
customers  and  their  weapon  systems  are  not  stationary:  customers  may  spend  a  certain 
time  in  close  proximity  to  one  fulfillment  center  and  then  move  closer  to  another 
fulfillment  center.  Future  work  may  model  mobile  customers  by  partitioning  the  network 
into  regions  and  implementing  a  region  index  that  specifies  where  particular  events 
occur.  Building  on  this,  we  recommend  that  future  work  formulates  business  rules  that 
specify  which  fulfillment  centers  provide  replacement  parts  to  which  customer  for  each 
time  period. 


3.  Further  Analysis 

In  terms  of  further  analysis,  we  foresee  great  utility  in  applying  our  formulation  to 
larger  multi-echelon  supply  networks.  As  it  stands,  our  analysis  only  considers  the 
notional  experiments  discussed  in  Chapter  III.  Real-world  weapon  systems  and  supply 
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networks  may  encompass  hundreds  of  echelons,  thousands  of  repairable  parts,  and 
millions  of  investment  dollars,  far  beyond  what  our  research  examines.  Future  work  may 
model  spare-part  management  of  an  actual  weapon  system  from  which  additional 
operational  insight  can  be  gained. 


49 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


50 


LIST  OF  REFERENCES 


[1]  C.  C.  Sherbrooke,  “Introduction,”  in  Optimal  Inventory  Modeling  of  Systems: 
Multi-Echelon  Techniques,  F.  S.  Hillier,  Ed.,  ed.  Boston,  MA:  Kluwer 
Academic  Publishers,  2004,  pp.  3-10. 

[2]  J.  A.  Muckstadt,  “A  continuous  time,  multi-echelon,  multi-item  system  with  time- 
based  service  level  constraints,”  m  Analysis  and  Algorithms  for  Service  Parts 
Supply  Chains,  T.  V.  Mikosch,  S.  M.  Robinson,  and  S.  I.  Resnick,  Eds.  New 
York:  Springer,  2005,  pp.  109-138. 

[3]  P.  Chandra,  “A  dynamic  distribution  model  with  warehouse  and  customer 
replenishment  requirements,”  J.  Oper.  Res.  Soc.,  vol.  44,  no.  7,  pp.  681-692,  Jul. 
1993. 

[4]  H.  E.  Lee,  “A  Multi-echelon  inventory  model  for  repairable  items  with  emergency 
lateral  transshipments,”  Management  Science,  vol.  33,  no.  10,  pp.  1302-1316, 

Oct.  1987. 

[5]  H.  Pirkul  and  V.  Jayaraman,  “Production,  transportation,  and  distribution 
planning  in  a  multi-commodity  tri-echelon  system,”  Transportation  Science,  vol. 
30,  no.  4,  pp.  291-302,  Nov.  1996. 

[6]  J.  Schneider  and  S.  Kirkpatrick,  “Overview  of  optimization  heuristics,”  in 
Stochastic  Optimization.  Leipzig,  Germany:  Springer  Berlin  Heidelberg,  2006, 
pp.  43-47. 

[7]  S.  Axsater,  “A  new  decision  rule  for  lateral  transshipments  in  inventory  systems,” 
Management  Science,  vol.  49,  no.  9,  pp.  1168-1179,  Sep.  2003. 

[8]  S.  C.  Graves,  “A  multi-echelon  inventory  model  for  a  repairable  item  with  one- 
for-one  replenishment,”  Management  Science,  vol.  31,  no.  10,  pp.  1247-1256, 
Oct.,  1985. 

[9]  P.  Tsiakis,  N.  Shah,  and  C.  C.  Pantelides,  “Design  of  multi-echelon  supply  chain 
networks  under  demand  uncertainty,”  Ind.  Eng.  Chem.  Res.,  vol.  40,  pp.  3585- 
3604,  Jul.  2001. 

[10]  K.  E.  Caggiano,  P.  L.  Jackson,  J.  A.  Muckstadt,  and  J.  A.  Rappold,  “Efficient 
computation  of  time-based  customer  service  levels  in  a  multi-item,  multi-echelon 
supply  chain:  A  practical  approach  for  inventory  optimization,”  Eur.  J.  Oper. 

Res.,  vol.  199,  pp.  744-749,  Dec.  2009. 


51 


[11]  R.  Ganeshan,  “Managing  supply  chain  inventories:  A  multiple  retailer,  one 
warehouse,  multiple  supplier  model,”  Int.  J.  Prod.  Econ.,  vol.  59,  pp.  341-354, 
1999. 

[12]  J.  J.  Neale  and  S.  P.  Willems,  “Managing  inventory  in  supply  chains  with 
nonstationary  demand,”  Interfaces,  vol.  39,  no.  5,  pp.  388-399,  Sep.-Oct.  2009. 

[13]  T.  lida,  “The  infinite  horizon  non-stationary  stochastic  multi-echelon  inventory 
problem  and  near-myopic  policies,”  Eur.  J.  Oper.  Res.,  vol.  134,  pp.  525-539, 
Nov.,  2001. 

[14]  D.  Simchi-Levi  and  Y.  Zhao,  “Three  generic  methods  for  evaluating  stochastic 
multi-echelon  inventory  systems,”  unpublished. 

[15]  M.  Ettl,  G.  E.  Eeigin,  G.  Y.  Ein,  and  D.  D.  Yao,  “A  supply  network  model  with 
base-stock  control  and  service  requirements,”  Oper.  Res.,  vol.  48,  no.  2,  pp.  216- 
232,  Mar.-Apr.,  2000. 

[16]  J.  Acimovic  and  S.  Graves,  “Making  better  fulfillment  decisions  on  the  fly  in  an 
online  retail  environment,”  unpublished. 

[17]  S.  Niranjan  and  F.  W.  Ciarallo,  “Base-stock  levels  in  multi-echelon  inventory 
systems  with  multiple  intermediary  product  demands,”  in  Industrial  Engineering 
Research  Conf,  Nashville,  TN,  2007,  pp.  298-303. 

[18]  GAMS.  2012.  General  Algebraic  Modeling  System.  [Online].  Available: 
http://www.gams.com/download 

[19]  GAMS.  2012.  CPEEX.  [Online].  Available: 
http://www.gams.com/dd/docs/solvers/cplex.pdf 


52 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Technical  Information  Center 
Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 
Naval  Postgraduate  School 
Monterey,  California 


53 


